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We discuss the role of many-body spin correlations in nonlinear optical response of a Fermi sea 
system with a deep impurity level. Due to the Hubbard repulsion between electrons at the impurity, 
the optical transitions between the impurity level and the Fermi sea states lead to an optically- 
induced Kondo effect. In particular, the third-order nonlinear optical susceptibility logarithmiclly 
diverges at the absorption threshold. The shape of the pump-probe spectrum is governed by the 
light-induced Kondo temperature, which can be tuned by varying the intensity and frequency of the 
O \ pump optical field. In the Kondo limit, corresponding to off- resonant pump excitation, the nonlinear 

absorption spectrum exhibits a narrow peak below the linear absorption onset. 



I. INTRODUCTION 



The role of many-body correlations in the nonlinear optical spectroscopy of Fermi sea systems has attracted much 
interest during the last decade [0-j3|. Due to the development of high quality ultrashort laser pulses, it has become 
possible to probe the elementary excitations of an interacting system on time scales shorter than the dephasing 
\ and relaxation times. In undoped semiconductors, it has been established that exciton-exciton interactions play a 
i dominant role in the coherent regime 0], Very recently, there has been a growing interest in studying the coherent 
| dynamics of Fermi sea (FS) systems at low temperatures 

In this paper, we discuss a new many-body effect in the nonlinear optical response of a FS system with a deep 
impurity level. Specifically, we will focus on nonlinear absorption due to optical transitions between a localized 
impurity level and the continuum of FS states. In the linear absorption, two prominent many-body effects has long 
been known in such systems First is the Mahan singularity due to the attractive interaction between the FS 

and the localized hole. Second is the Anderson orthogonality catastrophe due to the readjustment of the FS density 
profile during the optical transition. 

The nonlinear absorption comes from multiple transitions between impurity level and the Fermi sea. In particular, a 
number of different intermediate processes contribute to the third-order optical susceptibility x'' 3 '' IMI < What is crucial 
for us here is that, in the system under study, some of the intermediate states involve the doubly-occupied impurity 
• i-h , level. For example, the optical field can first cause a transition of a FS electron to the singly-occupied impurity level, 
1 which thus becomes doubly-occupied, and then excite both electrons from the impurity level to the conduction band. 
This is illustrated in Fig. 1(a). Important is that, while on the impurity, the two electrons experience a Hubbard 



repulsion. In fact, such a repulsion gives leads to a logarithmic divergence in ■ The origin of such an anomaly is 
intimately related to the Kondo effect fj^] . 

To be specific, we restrict ourselves to pump-probe spectroscopy, where a strong pump and a weak probe optical 
field are applied to the system and the optical polarization along the probe direction is measured. We only consider 
near-threshold absorption at zero temperature and assume that the pump frequency is tuned below the onset of 
optical transitions from the impurity level so that dephasing processes due to electron-electron and electron-phonon 
interactions are suppressed. 

We identify two distinct regimes which are characterized by the interplay between pump detuning and pump 
intensity. These regimes are somewhat analogous to the Kondo limit and mixed-valence regime in the "usual" Kondo 
systems. The shapes of nonlinear absorption spectra in these regimes are drastically different. The crossover between 
the two regimes is governed by a new energy scale - the lint-induced Kondo temperature, which can be tuned by 
varying the intensity and frequency of the pump. 

The paper is organized as follows. In Section II, we derive the Kondo-induced contribution to the third-order optical 
polarization that determines the absorption spectrum. In Section III, we obtain the absorption coefficient beyond 
v/ 3 ) using variational large N method. Section IV concludes the paper. 
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FIG. 1. Intermediate processes contributing to X ■ ( a ) Intermediate state with doubly-occupied impurity: excitation of a 
FS electron to the second impurity state followed by the excitation of the two impurity electrons to the FS. (b) Large U limit: 
Two transition channels are possible from states below the FS to the two spin-degenerate states of empty impurity, but only 
one channel from the singly-occupied impurity to states above the FS. 



II. KONDO ANOMALY IN X 



(3) 



We start with the Hamiltonian of the system: H tot — H + H\{t) + H 2 (t), where 



(1) 



ktr 



is the Hamiltonian in the absence of optical fields; here cj^ and d£ are conduction and localized electron creation 
operators, respectively, {h„ — S a d a ); Sk and Ed are the corresponding energies, and U is the Hubbard interaction (all 
energies are measured from the Fermi level). The coupling to the optical fields is described by the Hamiltonian 



Hi{t) = -Mi(t)ft + h.c, f t = 4<A> Mi(t) 



^ik-i-r—iujit 



(2) 



where T* is the optical transition operator. Here £i(t), k^, and u>i are the pump/probe electric field amplitude, 
direction and central frequency, respectively, and /i is the dipole matrix element (i = 1,2 denotes the probe and 
pump, respectively). The pump-probe polarization is obtained by expanding the optical polarization, ^(T 1 ), to the 
first order in H\ and keeping the terms propagating in the probe direction pT| : 



P(t) =iji i dt'M^t') 



(§{t)\fK{t,t')f^\§(t')) - (<f>(t')\f^lC(t',t)f\<S>(t)) 



(3) 



where K.(t,t') is the evolution operator for the Hamiltonian H + H2(t) and the state | <!>(£)) satisfies the Schrodinger 
equation 0t|$(t)) = [H + H 2 (t)]\$(t)) . 

The third order polarization is obtained by expanding K.(t,t') and up to the second order in H 2 . Below we 

consider sufficiently large values of U so that, in the absence of optical fields, the ground state of H, \Slo), represents 
a singly- occupied impurity and full FS. For large U, the doubly-occupied impurity states are energetically unfavorable 
and can be excluded from the expansion of the polarization (H) with respect to H 2 . The third-order pump- probe 



polarization then takes the form (t) 



iki -r— iu)\t 



with 



Qi(t, t') + Ql (*', t) + Q 2 (t, t') + Q 3 (t, t') 



(4) 
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where 

Ql(t,t') = - f dh [ dt 2 f(t 1 ,t 2 )F(t,t',t 1 ,t 2 ), 

J — OO J — OQ 

Q 2 (t,t') = -J / dt 2 J t dhfitt^WitMM'), 
Q 3 (t,t') = - f dh f dhf{ti,h)F{h,1?,t,h). 

J—OO J — OQ 

Here we denoted f(h,t 2 ) = ^(ii^^e 1 " 2 ^ 1 ^, and 

F(t,t',t u t 2 ) = (n \fe~ lHi - t ~ t ^f^e^ H< - t '- t ^fe- lH{tl - t ^f t \n ) 



A pqk 



a e - i ( £ p- e d)(t-t')-i(e k -e k , )(t'-ti )-i(sk -£d) *2) 



(5) 



(6) 



pqk'kAser'er 



^pqk'k = (^ol4 C P^ C qs d s d l' C k'<T' c L d ffl^O> = S\adsa>n a {l ~ n p ) [<Jpk<S qk '7Jg + <W<5pq<W(l ~ n fc)L 



(7) 



with n a = (flold^da-lfto) and rife = (Slo| c k<T c ko-|^o) (impurity occupation number is = n a — 1 here). Eqs. 
and (m) are valid in the large U limit corresponding to a singly-occupied impurity. For monochromatic optical fields, 
£i{t) = Si, the time integrals can be explicitly evaluated, 



Qi{t,t')+Q\{t,t') = ^^Te-^-^C-^l - n p ) 
pq 

Q 2 (t,t')= SiJ2 e ' i(e ''~ ed)it ~ t ' ) ^ - np) 
pq 

Q 3 (t,t')= -£ 2 2 ^(1 - n p ) 
pq 

yielding = P (3) + P^ 3) with 
pf 



27V n „ 



i(t-t')(l-n q ) 



(e p - e q )(e p - E d ) (e q -E d ) 

i{t-t')Nn q 2(1 -n q ) (1 - n^e 1 ^-^'-*') 

(e q - E d ) (e p - e q )(e p - E d ) (e q - E d ){e p - E d ) 

e i(e p -e q -u 2 )(t-t') Nfl(j + e -iu, 2 {t-t')(l _ n ^ 

(e q - E d ){e p - E d ) 



(1 - n p ) 



(8) 



(e p - e q )(s p - E d ) {e p - s d - uji)(s q - E d ) 



PP = {N-l)^£ 1 ElY, 



(1 - n p )n q 
£p — £d — Wi 



(s p - s q )(s p - E d ) {s p - s d - ui)(s q - E d ) 



(9) 
(10) 



where N is the impurity level degeneracy. Here we introduced the effective impurity level Ed = s d + ^2 ■ The first 
term, Pq , is the usual third-order polarization for spinless (N = 1) electrons pj| . The second term, P^, originates 
from the suppression, due to the Hubbard repulsion U, of the contributions from doubly-occupied impurity states. 
As indicated by the prefactor (N — 1), it comes from the additional intermediate states that are absent in the spinless 
case [see Fig 1(b)]. 

Consider the first term in Eq. (|l0|). The restriction of the sum over q to states below the Fermi level results in a 
logarithmic divergence in the absorption coefficient, a oc ImP, at the absorption threshold, u>i = —Ed'. 



ImP® = (N- l)po0(wi + £d)^r- In 

irouj 



2A 



D 



LO\ + Ed 



(11) 



where po = 7r£i/j 2 t/, 5uj = uj\ — lo 2 is the pump-probe detuning, and A = irg[j, 2 £ 2 is the energy width characterizing 
the pump intensity; D and g are the bandwidth and the density of states (per spin) at the Fermi level, respectively. 
Recalling that the linear absorption is determined by ImPW = p 9(uji + e d ), we see that it differs from Eq. (Jll|) by 



a factor ln| Ml ^ e | (setting for simplicity N — 2). In other words, ImP' 1 ) and IruP£? become comparable when 
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u>i + Ed = Slu + Ed r~j D exp 



ttSlo 



(12) 



We see that the perturbative expansion of the nonlinear optical polarization in terms of the optical fields breaks down 
even for weak pump intensities (i.e., small A). The above condition of its validity depends critically on the detuning 
of the pump frequency from the Fermi level. For off-resonant pump, such that the effective impurity level lies below 
the Fermi level, 123^1 = \ej\ ~ ^2 S> A, the relation (n3) can be written as 5oj + Ed ~ Tr- where 



T K = De* Ed/2A = Dexp 



\£d\ - u 2 



(13) 



This new energy scale can be associated with the Kondo temperature — an energy scale known to emerge from a 
spin-flip scattering of a FS electron by a magnetic impurity Remarkably, in our case, the Kondo temperature 
can be tuned by varying the frequency and intensity of the pump. In fact, the logarithmic divergence in Eq. ( |TT| ) is 
an indication of an optically-induced Kondo effect. 

Let us now turn to the second term in Eq. (|Io|). In fact, it represents the lowest order in the expansion of the linear 
polarization with impurity level shifted by Se — (N — l)/i 2 £f 



q Eg~E d 



(1-Tlp) 



Ed + Se - lui 



(14) 



The origin of Se can be understood by observing that, for monochromatic pump, the coupling between the FS and the 
impurity can be described by a time-independent Anderson Hamiltonian Ha with effective impurity level Ed — £d + ^2 
and hybridization parameter V = By virtue of this analogy, Se itself is the perturbative solution of the following 
equation for the self-energy part |0| 



E = E(Bb) = (N- \)^£l 



US ^(N-l£ln Ed - E °. 

: q -E d + E y n D ' 



(15) 



which determines the renormalization of the effective impurity energy, Ed-, to Ed = Ed — Eq [ |T3[ . Indeed, to the first 
order in the optical field, Eq. ( flit ) yields Eq = Se after omitting Eq in the rhs. 

The logarithmic divergence ( |il| ) indicates that near the absorption threshold, a nonperturbative treatment is nec- 
essary. Recall that, in the lowest order in vq, the attractive interaction vq between a localized hole and FS electrons 
also leads to a logarithmically diverging correction even in the linear absorption: SP^ 1 ' ~ PWjbj ln[D/(o;i + Ed)]- In 
the nonperturbative regime, 5pW ~ pW, this correction evolves into the Fermi edge singularity JTo| . The question 
is how the Kondo correction (|ll]) will evolve in the nonperturbative regime. This question is addressed in the next 
section. 



III. GENERAL SHAPE OF ABSORPTION SPECTRUM 



We start by discussing qualitatively the results and defer the details to the end of the section. 
It can be seen from the expression ( [l3| ) for Tjc that there is a well-defined critical pump intensity, 

A c = n 9f i 2 £l = |(| Ed |-u, 2 ). (16) 

Note that |e^| — ^2 is the pump detuning from the Fermi level. The shape of the nonlinear absorption spectrum 
depends sharply on the ratio between A and A c . The strong pump case, A > A c , is analogous to mixed-valence 
regime. In this regime, the Kondo correction ( |Tl| ) develops into a broad peak with width A and height pq. This is 
illustrated in Fig. 2(a). 

Much more delicate is the case A <C A c , which is analogous to the Kondo limit. The Kondo scale Tk is then much 
smaller than A, which is the case for well-below-resonance pump excitation, \ed\ - W2 > A. The impurity density 
of states in the Kondo limit is known |l3| to have two peaks well separated in energy by \Ed\ — |ed| — 0J2 ^> A (Ed 
is the effective level position). As a result, in the presence of the pump, the system sustains excitations originating 
from the beats between these peaks. These excitations can, in fact, assist the absorption of a probe photon. The 
corresponding condition for the probe frequency reads \E<i\ + u)\ ~ \^d\, or wi ~ u>2. Thus, in the Kondo limit, the 
absorption spectrum exhibits a narrow peak below the linear absorption onset. This is illustrated in Fig. 2(b). 
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FIG. 2. Schematic plot of the Kondo-absorption spectra vs probe frequency, (a) Mixed-valence regime: spectrum for strong 
pump intensity (thick line) compared with the perturbative result ( |ll| ) (dashed line) and the linear absorption spectrum (thin 
line), (b) Kondo limit: the nonlinear absorption spectrum exhibits a narrow peak below the linear absorption threshold. 

To calculate the shape of the below-thrcshold absorption peak, we adopt the large N variational wave-function 
method by following the approach of P4| . For monochromatic optical fields, the polarization (|3j) can be written as 



P = —fi£i 



G < (E -5lu) + G > (E + 6lu) 



(17) 



where 



G<(e) = (n\T^e-H A )- 1 T\n) 1 (18) 
and G > (e) is similar but with T < — > T>. In the leading order in TV -1 , the ground state is given by 

\n)=A(\Q)+J2n q a q \q,l)), (19) 



where |q, 1) = iV" 1 / 2 4cq<r|0) (|0) stands for the full FS). The coefficients A and a*; are found by minimizing Ha 
in this basis, 



A 2 — 1 - rid, a q 



Ed — e« — Eq' 



where 



rid = 1 



7TE d 

NA 



(20) 



(21) 



is the impurity occupation |lj,|l3(] (NA is finite in the large N limit). Using that T|f2) = vNj^^ | q, 1) , the Green 
function G < (e) takes the form 



G<(e) = NA 2 ^n p n q (v>,l 



e-H A 



q,l , 



(22) 



5 



where the matrix element in the rhs is obtained, to the leading order in TV , by inverting Ha in the above basis set, 



P,l 



q,l 



with E(e) given by Eq. (|15|). We then obtain 



G<(E -5u) = - 



E d 



H{E - 5lu) 



S(e) (e + e 9 -£ d )(e 



|s(g -MI 2 

E -5w- £(£ - <M 



(23) 



(24) 



Since T,(Eq) = Eq [see Eq. (|l5|)], the second term has a pole at 8oj = which gives rise to a resonance. The N 1 
correction gives a finite resonance width A. The first term gives nonresonant contribution to absorption for 5lu < Ed- 
In a similar way, it can be shown that contribution from G > (Eo + Sui) is suppressed by factor iV^ 1 and is nonresonant. 
Near the resonance, using that [dT,(Eo) / OEq — = rid — 1 PI, we obtain 



lmP K = po 



El{l~n d f 



Slo 2 + A 2 

For A -C A c , corresponding to the Kondo limit, we have 1 — n d — ttTk/NA and Eq ~ Ed, so that 



ImPx ~ po 



1 



( nE d T K Y 

V NA J Suj 2 + A 2 ' 



(25) 



(26) 



In this case, ( p6[ ) describes the narrow below-threshold peak [see Fig. 2(b)]. The factor (1 — rid) 2 has the physical 
meaning of a product of populations of electrons in the narrow peak of the impurity spectral function (Kondo 
resonance) and "holes" in the wide peak (centered at Ed below the Fermi level). In the opposite case, A > A c , 
corresponding to mixed-valence regime, we have 1 — rid ~ 1 and Eq ~ NA. Then the polarization 



lmP K = Po 



N 2 A 2 
Suj 2 + A 2 



(27) 



describes the absorption peak in Fig. 2(a). 



IV. COCLUSION 



Note that, although we considered here, for simplicity, the limit of singly occupied impurity level in the ground 
state, the Kondo-absorption can take place even if the impurity is doubly occupied. Indeed, after the probe excites an 
impurity electron, the spin-flip scattering of FS electrons with the remaining impurity electron will lead to the Kondo 
resonance in the final state of the transition. In this case, however, the Kondo effect should show up in the fifth-order 
polarization. 

A feasible system in which the proposed effect might be observed is, e.g., GaAs/AlGaAs superlattice delta-doped 
with Si donors located in the barrier. The role of impurity in this system is played by a shallow acceptor, e.g., Be. 
MBE growth technology allows one to vary the quantum well width and to place acceptors right in the middle of each 
quantum well p5[ . In quantum wells, the valence band is only doubly degenerate with respect to the total angular 
momentum J. Thus, such a system emulates the large U limit considered here. The dipole matrix element for acceptor 
to conduction band transitions can be estimated as p. ~ Paa, where fj,Q is the interband matrix element and a is the 
size of the acceptor wave function. For typical excitation intensities g), the parameter A ranges on the meV scale 
resulting in Tjc ~ A for the pump detuning of several meV. 

Finally, let us discuss the effect of a finite duration of the pump pulse, t. Our result for remains unchanged 
if t is longer than K/Tk- If r < K/Tk, then r will serve as a cutoff of the logarithmic divergence in (pM|), and the 

Kondo correction will depend on the parameters of the pump £2 and t as follows: ImP^- oc £| \t\(Dt / K) . In the 
non-perturbative regime, our basic assumption was that, for monochromatic pump, the system maps onto the ground 
state of the Anderson Hamiltonian. Our results apply if the pump is turned on slowly on a time scale longer than 
Ti/Tk- For shorter pulse duration, the build up of the optically-induced Kondo effect will be determined by the 
dynamics of FS excitations. 
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